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Abstract 

The reconstruction of reticulate evolutionary histories in plants is still a major methodological challenge. Sequences of the 
ITS nrDNA are a popular marker to analyze hybrid relationships, but variation of this multicopy spacer region is affected by 
concerted evolution, high intraindividual polymorphism, and shifts in mode of reproduction. The relevance of changes in 
secondary structure is still under dispute. We aim to shed light on the extent of polymorphism within and between sexual 
species and their putative natural as well as synthetic hybrid derivatives in the Ranunculus auricomus complex to test 
morphology-based hypotheses of hybrid origin and parentage of taxa. We employed direct sequencing of ITS nrDNA from 
68 individuals representing three sexuals, their synthetic hybrids and one sympatric natural apomict, as well as cloning of 
ITS copies in four representative individuals, RNA secondary structure analysis, and landmark geometric morphometric 
analysis on leaves. Phylogenetic network analyses indicate additivity of parental ITS variants in both synthetic and natural 
hybrids. The triploid synthetic hybrids are genetically much closer to their maternal progenitors, probably due to ploidy 
dosage effects, although exhibiting a paternal-like leaf morphology. The natural hybrids are genetically and morphologically 
closer to the putative paternal progenitor species. Secondary structures of ITS1-5.8S-ITS2 were rather conserved in all taxa. 
The observed similarities in ITS polymorphisms suggest that the natural apomict R. variabilis is an ancient hybrid of the 
diploid sexual species R. notabilis and the sexual species R. cassubicifolius. The additivity pattern shared by R. variabilis and 
the synthetic hybrids supports an evolutionary and biogeographical scenario that R. variabilis originated from ancient 
hybridization. Concerted evolution of ITS copies in R. variabilis is incomplete, probably due to a shift to asexual 
reproduction. Under the condition of comprehensive inter- and intraspecific sampling, ITS polymorphisms are powerful for 
elucidating reticulate evolutionary histories. 
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introduction 

Hybridization and polyploidy are creative evolutionary forces in 
plant radiations, and most or perhaps all angiosperms are either 
polyploid or of ancient polyploid origin [1]. High morphological 
plasticity and large ecological amplitude of allopolyploids could be 
derived from an increased propensity for epigenetic variation 
following interspecific hybridization [2]. Rieseberg and collabora- 
tors [3] proposed that ancestral phenotypic and genomic traits of 
recent hybrids might be recognized by comparison with experi- 
mental hybrids. Hybrid speciation is facihtated in plants through 
subsequent polyploidization [4] and can be connected to a 
transition from sexuality to apomixis, i.e., asexual reproduction via 
seed [5] . Nonetheless, the reconstruction of hybrid relationships 
and parentage in natural systems is stiU a methodological 
challenge. 

The internal transcribed spacer (ITS) of the nrDNA has been 
established as a standard molecular marker to infer generic and 
interspecific relationships in flowering plants [6-8] and to infer 



hybridization events in other eukaryotes including algae, hetero- 
trophic protists, invertebrates and vertebrates [9-15]. ITS nrDNA 
has proved to be a helpful non-coding marker to infer 
hybridization events because of its biparental inheritance and its 
occurrence in hundreds to thousands of copies within a single 
genome [8,16,17]. The majority of systematics studies aim at 
reconstruction of phylogenetic relationships using ITS in combi- 
nation with maternally inherited plastid marker (e.g., [18-23]). 

According to Song and collaborators [24], a genome of each 
plant individual contains 35 ITS variants on average. The 
homogenization of different progenitor ITS variants in sexual 
hybrids through concerted evolution is a well-documented 
phenomenon [25—27]. In apomictic plant lineages, the homoge- 
nizing effects of meiosis are lacking, and thus heterozygous ITS 
polymorphisms inherited from hybrid ancestry may persist over 
several generations [21,28-30]. Only few studies have analyzed 
ITS polymorphisms in apomictic hybrids [21,28,29,31,32]. The 
features of secondary structure of ITS in apomictic plants, 
however, are still unknown. 
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Figure 1. Distribution of the analyzed populations of the Ranunculus auricomus comp\sx. Dots represent sampled populations, the areas 
with the same color show the entire distribution of the respective taxon within Central Europe (after [49]). 
doi:10.1371/journal.pone.0103003.g001 



Apart from the lasting popularity of this marker, there is stiU a 
lively debate whether a direct link exists between ITS structural 
features like compensatory base changes (CBCs) and Ernst Mayr's 
biological species concept [33-35] . CBCs are mutations that occur 
in a primary RNA transcript, whereby both nucleotides paired in 
the secondary configuration of the ITS transcript mutate so that 
their bond is retained (e.g., G-C mutates to A-U). A hemi-CBC 
(hCBC) is the mutation of one of the two nucleotides while 
maintaining the nucleotide bond. Non-compensatory base changes 
alter the secondary structure by disabling pairing within any of the 
helical regions (i.e., within the internal paired regions). According 
to [34], the occurrence of CBCs in ITS2 secondary structure 
between two sexual species is correlated to unsuccessful interspe- 
cific mating. This assumption is based solely on coincidence, i.e., 
similar rates of evolutionary changes affecting both the ITS 
sequences and the genes involved in sexual reproduction [33,36]. 
Hence, the appearance of a CBC between species could 
theoretically serve as a criterion for delimitation of sexual species. 
Reflecting the difficulty associated with multiple copies of the ITS 
within a genome, analysis of an extensive ITS2 dataset obtained 
by high-throughput sequencing in plants suggests that intrage- 
nomic CBCs are unlikely. Thus, CBCs can potentially serve as 
species-specific markers [24,34]. 

Despite the wide usage of the ITS region, only a few studies 
have discussed obvious difficulties regarding the presence of both 
functional and putative non-functional copies (e.g., [37^1]). The 
homogenization of rRNA genes through gene conversion prevents 
these from accumulating mutations [42]. The most frequently 
employed criteria to recognize putative non-functional ITS copies 
(also referred to as "pseudogenes") are nucleotide substitutions in 
highly conserved motifs (especially within 5.8S sequence regions), 
occurrence of methylation-induced substitutions (C ^ T, G ^ A), 



decreased or varying GC content, incorrect secondary structure 
folding and incongruent phylogenetic placement of conspecific 
ribotypes [21,37]. For this reason, frequent occurrence of putative 
non-functional copies may be a reliable hint for the incompleteness 
or even suppression of concerted evolution of the ITS1-5.8S-ITS2 
copies [21,39,43,44]. 

As in many other plant genera [45], hybridization and 
polyploidization have shaped the evolutionary history of butter- 
cups {Ranunculus) [29,46-51]. The Euro-Siberian R. auricomus 
complex, which belongs to one clade classified as R. sect. 
Auricomus [51], comprises about 800 polyploid apomictic 
microspecies and only four sexual species [52]. Here we investigate 
three sexual species {R. carpaticola So6, R. cassuhicifolius W. 
Koch, R. notahilis HorancU and Guterm.) and one apomictic 
species (R. vaiiahilis Horandl and Guterm.) from Central Europe. 
These species are taxonomicaUy subdivided into two morpholog- 
ical collective groups: the auricomus group { = R. notahilis, R. 
variabilis) and the cassubicus group {~R. carpaticola, R. 
cassuhicifolius) [53]. According to population genetic cUstances, 
these two groups separated about 0.9 Ma ago [54]. Autotetraploid 
R. cassuhicifolius and diploid R. carpaticola are woodland species 
of the pre-Alps and Carpathians, respectively, whereas diploid R. 
notahilis is endemic in a small lowland area in southeastern Austria 
and occurs in woods and meadows; allotetraploid R. variahilis is 
widespread in lowlands of Austria, sympatric with R. notahilis and 
R. cassuhicifolius, but colonizing mostly anthropogenic meadows 
(Figure 1). AUozyme studies and morphometries based on leaves 
and fruit characters suggested that diploid R. notahilis (auricomus 
lowland group) is a putative progenitor of tetraploid R. variahilis 
[55]. Surprisingly, experimental crosses between 4x R. cassuhici- 
folius and 2x R. notahilis, employing the latter as pollen donor, 
resulted in triploid synthetic hybrids that resembled R. variahilis in 
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Table 1. List of directly sequenced individuals. 



taxon 


herbarium collection 
No. or garden 
accession No. 


population 


location/GPS coordinates 


ploidy 
level 


sexual or 
apomictic 


GenBank 
accession No. 


R. notabilis 


7223-6/-10 


N0T2 


AT, Strem/47° 2' 58.0" N; 
Id ZD J./ b 


2x 


sex''' 


KF671993-KF671994 


R. notabilis 


561 3-1 /-2 


N0T3 


AT, Strem/47=^ 2' 58.0" N; 
16^ 26' 3.7" E 


2x 


sex^' 


KF671996-KF671997 


R. notabilis 


7224-15 


N0T3 


AT, Strem/47 0' 56.0" N; 
Id z6 Jo.d b 


2x 


sex^' 


KF671995 


R. notabilis 


9609-2/-5/-9 


N0T3 


AT, Strem/47=^ 2' 58.0" N; 
16" 26' 3.7" E 


2x 


sex^' 


KF671999-KF672001 


R. notabilis 


5615 


N0T4 


AT, Strem/47° 2' 58.0" N; 
16"^ 26' 3.7" E 


2x 


sex^' 


KF671998 


R. variabilis 


7221-13/-6/-9 


BURG2 


AT, Strem/47 2' 51.3" N; 
Id 25 DD.U b 


4x 


apo^' 


KF672008-KF672010 


R. variabilis 


8210-5/-7/-22/-23/-25 


BUCKl 


AT, Spratzau/47" 35' 
iy.5 N, 16 12 li.o b 


4x 


apo"' 


KF672002-KF672006 


R. variabilis 


8212-23 


BUCK2 


AT, Spratzau/47 33' 
25. J N, 16 18 45.2 b 


4x 


apo^' 


KF672007 


R. variabilis 


821 4-1 /-2 


WALDl 


AT, Rosenburg/48° 37' 

52.0 N; 1 5 38 /.3 b 


4x 


apo"' 


KF672017-KF672018 


R. variabilis 


8215-13 


WALD2 


AT, Wegscheid/48= 36' 
25.1 N; 1 5 28 59.9 b 


4x 


apo^' 


KF672011 


R. variabilis 


8216-26/-33 


WALD3 


AT, Zwettl/48" 35' 
42./ N, 1 5 10 2.6 b 


4x 


apo"' 


KF672015-KF672016 


R. variabilis 


8218-24 


MUHL2 


AT, Kefermarkt/48' 

26 31.6 N; 14 33 6.8 b 


4x 


apo"' 


KF672012 


R. variabilis 


8222-22/-25 


MUHL3 


AT, Reichenau/48" 27' 
25./ N; 14 21 24.2 b 


4x 


apo^' 


KF672013-KF672014 


R. cassubicifolius 


8472-1 /-3/-10/-1 9 


YBBl 


AT, Schwarzois/47° 54' 
48.0 N; 14 56 3/.4 b 


4x 


sex"' 


KF672026-KF672029 


R. cassubicifolius 


8473-24/-25 


YBB2 


AT, Schwarzoi5/47" 52' 

29.5 N; 14 55 59./ b 


4x 


sex"' 


KF672030-KF672031 


R. cassubicifolius 


A7, 17, 25, 35 


8473 X8473 


DE, Gottingen, 
experimental crossing 


4x 


sex*^' 


KF672032-KF672035 


R. carpaticola 


8483-1 /-4/-5/-6/-10 


REVl 


SK, Revuca/48= 41' 
iU.9 N; 20 / 34.9 b 


2x 


sex"' 


KF672019-KF672025 


R. carpaticola 


8486-27-4 


REV2 


SK, Revuca/48 41' 
21 .9" N; 20' 5' 45.4" E 


2x 


sex"' 


KF672023-KF672024 


R. carpaticola x notabilis 


F2, 4, 7, 1 0, 90 


8483x7224 


DE, Gottingen, 
experimental crossing 


2x 


sex" 


KF671968-KF671972 


R. carpaticola x notabilis 


JIO, 20, 22, 25, 33, 35 


8483x7224 


DE, Gottingen, 
experimental crossing 


2x 


sex^' 


KF671987-KF671992 


R. cassubicifolius x notabilis 


Gl, 2, 5, 7, 9, 11, 15-20 


8472x7220 


DE, Gottingen, 
experimental crossing 


3x 


sex"' 


KF671973-KF671984 


R. cassubicifolius x notabilis 


H4, 5 


8473x7223 


DE, Gottingen, 
experimental crossing 


3x 


sex* 


KF671985-KF671986 



"[96]. 
"'[46]. 

^'functional sexual seed, but with aposporous initials {for details see [58]). 
'^'functional sexual seed, but with low rates of aposporous seeds (for details see [58]). 
doi:l 0.1 371 /journal.pone.Ol 03003.t001 



basal leaf shape. Since the geographical range of R. variahilis 
spans over that of the sexual species (Figure 1), we hypothesize an 
allopolyploid, Pleistocene origin of i?. variabilis from hybrids of i?. 
cassubicifolius and R. notabilis. The fourth known sexual species 
of the R. auricomus complex (i.e., R. marsicus Guss. et Ten.) was 
not taken into account because its distribution area is restricted to 
the Central Apennines in Italy, which is geographically isolated 



from all other species under study [56]. Moreover, it is tetraploid 
and hexaploid and therefore must have another evolutionary 
origin. All sexual species of the complex belong to one clade but 
the internal relationships of this clade are not yet well resolved 
[57]. However, the other species of this clade (i.e., of section 
Auricomus) occur in other continents, and therefore interspecific 
hybridization is unlikely. 
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Table 2. List of cloned individuals. 





specimen 

(= clone library) 


clone 
identifiers 


population 


location/GPS 
coordinates 


ploidy 
level 


sexual or 
apomictic 


GenBank 
accession No. 


R. notabilis 7223-10 


N01-N20 


N0T2 


AT, Strem/47° 2' 58.0" N; 
16= 26' 3.7" E 


2x 


sex"' 


KF672076-KF672095 


ft. variabilis 8210-5 


vol -VI 9 


BUCK1 


AT, Spratzau/47° 35' 
19.5" N; 16° 12' 13.8" E 


4x 


apo" 


KF672096-KF672114 


ft. cassubicifolius 8472-1 


S01-S20 


YBB1 


AT, Schwarzois/47° 

54' 48.0" N; 14= 56' 37.4" E 


4x 


sex"' 


KF672036-KF672055 


ft. cassubicifolius x notabilis Gil 


X01-X20 


8472x7220 


DE, Gottingen, 
experimental crossing 


3x 


sex"' 


KF672056-KF672075 



"[96]. 
"'[46]. 

"'functional sexual seed, but with low rates of aposporous seeds (for details see [58]). 
doi:1 0.1 371 /journal.pone.Ol 03003.t002 



Here we aim to shed light on intraspecific ITS nrDNA 
variability in closely related, but morphologically diversified taxa 
of the Ranunculus auricomus complex in Central Europe. 
Therefore, an extended sampling of ITS ribotypes on the 
intraspecific level will enable us to recover additivity patterns in 
hybrid genomes to better understand complicated reticulate 
relationships among sexuals and apomicts. We want to address 
the following specific questions: (1) Do ITS polymorphisms 
characterize the sexual species and the hybrids? (2) Are ITS 
polymorphisms observed in artificial hybrids similar to those in 
natural tetraploid apomictic R. variabilis, indicating allopolyploid 
origin? (3) Are the ITS data congruent with morphometric data in 
the putative parent-hybrid relationships? (4) Is secondary structure 
of the entire ITS1-5.8S-ITS2 region informative about species 
differentiation, concerted evolution and formation of putative non- 
functional copies? 

Materials and Methods 

Ethics statement 

The collections used for this study did not involve endangered 
or protected species and no specific permissions were required for 
sampling activities in these locations. 

Plant material 

DNA samples were collected from R. carpaticola, R. cassubi- 
cifolius, R. notabilis, and R. variabilis (Table 1). Additionally, 
homoploid and heteroploid experimental hybrids (Fi offspring) 
obtained from hand-made crosses oiR. carpaticola x notabilis (1 1 
individuals), as well as R. cassubicifolius x notabilis (14 
individuals) were studied. In all crosses, R. notabilis was used as 
pollen donor; homoploid and heteroploid crosses were performed 
to reconstruct pathways to apomixis [58]. Ploidy levels of 
experimentally produced hybrids were assessed in mitotic cell 
divisions on root tips treated with a-bromonaphthalene during 
6 h, 3 h at RT and 3 h at 5°C, and fixed 24-48 h in 3:1 
(ethanohacetic acid). Then, root tips were hydrolyzed 10 min in 
IN HCl at 60°C and stained with Schiffs reagent (C20H21N3SO3). 
On average, five metaphase cells per plant were examined with a 
Leica B5500 microscope (Leica, Wetzlar, Germany). For prove- 
nance, accession codes, ploidy and method of reproduction of the 
sampled taxa see Table 1 and Figure 1 . The DNA sequencing was 
performed on 68 individuals. Plants were cultivated at the 
Botanical Garden of the University of Gottingen, Germany. 
Herbarium material was provided by the herbarium of the 



University of Vienna, Austria (WU). Four individuals were selected 
for cloning and 19-20 clones per individual were sequenced, 79 in 
total. A triploid individual (Gil) was selected rather than a diploid 
one because the maternal species from this cross occurs in the 
same geographical region as the putative natural hybrid (Figure 1). 
Cloned R. cassubicifolius 8472-1 represents the population used 
for crossing experiments. Another cloned individual R. notabilis 
7223-10 was taken from an isolated population to minimize 
presumed secondary introgressions from R. variabilis [55]. 
Additionally, we cloned the individual R. variabilis 8210-5 
belonging to a wild population geographically localized between 
those sexual species described above. 

DNA extraction, PGR amplification, cloning and 
sequencing 

Genomic DNA was extracted from sUica-dried, herbarium or 
frozen (— 80°C) plant material (stem leaves) using Tissue-Lyser II 
Qiagen and Invisorb Spin Plant Minikit (Invitek) following the 
manufacturer's instructions. PGR of the ITS1-5.8S-ITS2 region 
was performed according to [59] using primers 18sF and 26sR 
(after [60]); the PGR amplicons were cleaned with the Invisorb 
Spin PCRapid Kit (Invitek, Berlin, Germany). PGR products of a 
subset of four individuals {Ranunculus notabilis 7223-10, R. 
cassubicifolius 8472-1, R. cassubicifolius x notabilis Gil, R. 
variabilis 8210-5; Table 2) were cloned using the TOPO TA 
cloning kit (Invitrogen, Garlsbad, GA, USA) and the pGR2.1- 
TOPO vector. Ligation products were integrated into competent 
cells of Escherichia coli TOP 10, as supplied by the manufacturer. 
Plasmid DNA was extracted and purified with a NucleoSpin- 
Plasmid kit (Macherey and Nagel, Dtiren, Germany) following the 
manufacturer's instructions. Sequencing reactions were performed 
with a Dye Terminator Gycle Sequencing v3.1 kit (Applied 
Biosystems, Darmstadt, Germany) and an ABI Prism 3100 
(Applied Biosystems) automated sequencer. PGR primers were 
used for sequencing. Sequences were processed using the sequence 
analysis program Geneious Pro version 5.6.4 [61] and deposited in 
GenBank (KF671968- KF672114). 

Sequence alignment and data analysis 

Sequences comprising ITS1-5.8S-ITS2 were aligned with 
Glustal W [62] and further adjusted manually in BioEdit 7.0.9.0 
[63]. The online ITS2 database (http://its2.bioapps.biozentrum. 
uni-wuerzburg.de/) was used for annotation of the ITS2 region, 
whereas the ITSl and 5.8S were annotated by comparison to 
published secondary structures (e.g., [64,65]). The 5.8S region of 
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interspecific crossing scheme: 

F (2x) = Ranunculus carpaticol a (2x) x notabilis (2x) 
J (2x) = Ranunculus carpaticol a (2x) x notabilis (2x) 

G (3x) = Ranunculus cassubicifoliu s (4x) x notabilis (2x) 
H (3x) = Ranunculus cassubicifoliu s (4x) x notabilis (2x) 



auricomus morphotype: ^^i^ 
cassubicus morphotype: 



R. notabilis - "N" 

R. variabilis - "V" 

. uaticola - "P" 

R. cas subicifolius - "S' 



Figure 2. NeighborNet analysis of interspecific ITS1+ITS2 variability within the Ranunculus auricomus comp\ex. NeighborNet analysis of 
all ITS1+ITS2 sequences obtained by direct sequencing of the studied individuals. The spring leaf silhouettes illustrate the main phenotypic 
differences between the two morphotypes: the auricomus-morphotype (characteristic for R. notabilis, R. variabilis and the synthetic hybrid R. 
cassubicifolius x notabilis) and the cassub/cus-morphotype (i.e., R. carpaticola and R. cassubicifolius). Individuals belonging to R. carpaticola are 
marked as "P", R. cassubicifolius as "S", R. notabilis as "N" and R. variabilis as "V", respectively. Identical sequences representing the same ribotype are 
listed in Figure 3. Bootstrap values are given for the main clusters. 
doi:10.1371/journal.pone.0103003.g002 



directly sequenced accessions was completely conserved across the 
entire dataset, hence we excluded it from subsequent analyses and 
the final alignment consisted of ITS1+ITS2 regions (501 bp in 
total). The polymorphic sites were manually determined in BioEdit 
7.0.9.0 [63] and visualized separately for ITSl and ITS2 using an 
online apphcation to create DNA weblogos (http://weblogo. 



berkeley.edu/). Polymorphic positions of ITSl, 5.8S and ITS2 
were annotated using a single secondary structure model of R. 
notabilis (voucher specimen 5613-1) as uniform template (Figure 
SI). In order to visualize the reticulate relationships among the 
studied species, NeighborNet analysis was employed in Split- 
sTree4 4.10 [66], applying uncorrected P distances and ambigu- 
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ITS1 



ITS2 



[ploidy 1 
reprod. mode I 



cDooojcNor^mco-^r^ 



^cDOOT-o^oomo 

OOOT-T-OviCNjm^CD 



8483-AS1/-AS5 


T C ACCCC 


C 


cm 


G 


C 


CGG 


GC 


8483-AS4 r 2x 1 
8483-1 0/-AS6 sexual 


YC ACCCC 


C 


CT 


G 


C 


CGG 


GC 


YC ACCCC 


C 


C T 


G 


c 


CGG 


GC 


8486-2 *■ 


YC ACCCC 


C 


cM 


G 


c 


CGG 


GC 


8486-4 


r c d r r r c 


r 




c 
\j 


r 


r G c 




8472-1 /-3/-1 9 


T C ACCCC 


C 


cl 


G 


c 


CGG 


GC 


8472-10 


T C ACCCC 


C 


cl 


G 


c 


CGG 


GC 


8473-24 r 4x "1 


YC ACCCC 


c 


cl 


G 


c 


CGG 


GC 


8473-25 [sexualj 


T C ACCCC 


c 


c| 


G 


c 


CGG 


GC 


A7, A25, A35 


T C ACCCC 


c 


C T 


G 


c 


CGG 


GC 


A17 


T C ACCCC 


c 


C T 


G 


c 


CGG 


GC 


5613-1 


C C A T T T C 


T 


cM 


G 


c 


CGG 


GC 


5613-2 


C C A T T TMT 


C T 


G 


c 


CGG 


GC 


5615 


C CWT T T C T 


C T 


G 


YCGG 


GC 


7223-6 


CCATTTCT 


C T 


G 


C C G G 


GC 


7223- 10 r ^'^ 1 

1 sexual 1 

7224- 15 L J 


CCATTTCT 
CCATTTCT 


yI 


G 
G 


YMK KGC 

CCGGGC 


9609-2 


C C A T T TMT 


cE 


Q 


CCGG 


GC 


9609-5 


C C A T T T C 


T 


C T 


\J 


YMK KGC 


9609-9 


C C A T T T C 


T 


cm 




YA T T 


G Y 


F2A, F4A, F7A, F10A, F90 


YCA YY YC 


Y 


C T 


G 


C 


CGG 


11 


J10A, J20A, J22A [ 2x 1 


YCA YY YC 


YC T 


G 


C 


CGG 




J25A, J33A IsexuaJ 


Y YA Y Y YC 


YCEJ 




C 


CGG 


GC 


J35A 


YCAYY YC 


Y 


cm 


G 


C 


CGG 


GC 


Gil, G16A, G19A, G7A 


YCA Y Y YC 


Y 


cES 


G 


C 


CGG 




G20, G17A, G5, G9A, G18A 


YC A Y Y YC 


YCg 


G 


C 


CGG 


Vi 


r 3x 1 


YC A Y Y YC 


YC T 


G 
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polymorphic sites shared by R. notabilis + R. variabilis 
(R. notabilis + R. variabilis + synthetic hybrids = A) 
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Figure 3. Polymorphisms of ITS1+ITS2 witKiin tKie Ranunculus aun'comus comp\ex as detected by direct sequencing. A summary of all 
ITS1+ITS2 polymorphic sites (red and white letters) detected through direct sequencing. The three polymorphic sites which are shared by R. notabilis, 
its synthetic hybrid and the putative natural descendent (i.e., R. variabilis) are highlighted by white letters on black background and underlaid by 
green stripes. Those shared just by R notabilis and R. variabilis are underlaid by grey stripes. 
doi:l 0.1 371/journal.pone.01 03003.g003 



ities handled as average. Bootstrap support values for internal 
splits were calculated with 1000 replicates. Minimum energy 
modcLs of the secondar)' structures for the ITS1-5.8S-ITS2 rRNA 
transcripts were obtained using RNAstructure 5.3 [67]. For 
comparison, the 5.8S secondary structures from all cloned ITS 
nrDNA sequences were homologicaUy modeled based on the R. 
notabilis 5613-1 template using the online ITS2 database (http:// 
its2.bioapps.biozentrum.uni-wuerzburg.de/) and subsequently 
aligned in 4SALE 1.7 [68]. The final graphical output was 
generated by Varna 3.8 [69] and edited in CorelDra^\- Grapliics 
Suite X3. The GC contents within ITS1-5.8S-ITS2 regions were 
computed using MEGA 5.1. [70] and the resulting percentages for 
each clone were plotted as a 3D scatter plot in PAST 2.17c [71]. 
The cloned sequences were checked for putative non-functional 
ITS copies, chimeras and recombinants. First, the presence of 
angiosperm conserved motifs was determined as proposed by [40] 
in BioEdit 7.0.9.0 [63]. We focused especially on the nucleotide 
polymorphisms in the 5.8S region. This sequence region is stricdy 
conserved among all direcdy sequenced individuals and the single 
nucleotide polymorphisms observed in clones may be a hint for 
putative PGR artifacts or non-functional ITS copies. The 
alignment containing ITS1-ITS2 was checked for recombinants 
and chimeras using the program RDP 4.17 [72]. No chimeric 
sequences or recombinants were detected. The DnaSP 5.10 
software [73] was then used to delimitate ITS1-ITS2 ribotypes. 
Ribotype networks were constructed by statistical parsimony (with 
a parsimony probability set to 95%) in TGS 1.21 [74]. Finally, 
average Kimura 2-parameter genetic distances within and among 
groups were computed in MEGA 5.1. [70]. 

Landmark geometric morphometries 

In order to quantify the phenotypic variation and test for 
possible overlaps among the four genetically different taxa of the 
fflMncomiM-morphotype {R. notabilis, R. carpalicola x notabilis, R. 
cassubicifolius x notabilis, R. variabilis) we employed 2D- 
landmark geometric morphometric analyses. The auricomus 
morphotype shows a heterophyllous sequence of basal leaves with 
deeply divided leaf blades formed during the flowering period 
(basal leaves 3, 4, and 5 sensu [75], whereas in the cassuhicus 
morphotype the basal leaves formed during and after the flowering 
period are undivided (basal leaves 2 and 3 sensu [76]). In both 
groups the above-mentioned basal leaves of the flowering period 
are referred to as spring leaves. These taxonomicaUy most 
informative spring basal leaves (comp. [53]) were digitized from 
the above individuals employed for molecular investigation (or 
their respective populations) using a book2net Kiosk scanner. 
Subsequendy, 2D landmark data were obtained from leaf outlines. 
In total, data from 95 individuals were obtained for morphometric 
analyses (28 leaves of R. notabilis, 26 of the synthetic hybrids and 
41 oi R. variabilis). Fourteen fixed 2D landmarks were depicted 
on the leaf outline using TpsDig2 2.16 [77]. The landmarks were 
localized on the tips of major lobes and in the major incisions. The 
program TpsRelw 1.49 [78] was used to standardize the 
morphometric data via generahzed Procrustes analysis (GPA) 
[79]. Prior to subsequent multivariate analyses, landmark config- 
urations were bilaterally symmetrized according to [80] in the 
program PAST 2.17c [71]. The principal component analysis 
(PGA) of the geometric morphometric data (i.e., Relative Warps 



Analysis, RWA [79]) was conducted on the entire set of 95 leaves 
in TpsRelw 1.49 [78]. Herein we also reconstructed the 
hypothetical shape of mean landmark configuration within each 
group. Phenetic distances among the three compared taxa were 
expressed as Procrustes distances among the group mean 
landmark configurations (as computed in TpsSmaU 1.20 [81]). 
The statistical significance of difierences between groups was 
computed by two-group permutation tests (9999 permutations, 
Euclidean metric) using PAST 2.17c [71], based on specimen 
scores on the first 12 relative warps (accounting for 99.99% of the 
total variability in shape data). 

Results 

Inter- and intraspecific variability of ITS1-ITS2 inferred 
from direct sequencing 

NeighborNet analysis of ITS sequences (Figure 2) indicated a 
clear genetic separation of representatives of both morphotypes, 
the auricomus (dissected basal leaves, R. notabilis, R. variabilis) 
and the cassuhicus (undivided basal leaves, R. carpaticola, R. 
cassubicifolius). Both clusters were well supported but differed in 
their internal variability. The cassuhicus cluster appeared to be 
more homogeneous, i.e., diploid R. carpaticola and tetraploid R. 
cassubicifolius were rather intermixed than well separated from 
each other. The genetic similarity between the auricomus species 
was slightly lower (0.995) than between the cassuhicus species 
(0.997), and the mean genetic similarity between the both clades 
was even lower (0.982). The individuals within the auricomus clade 
exhibited a higher variability, consisting of several splits containing 
R. notabilis (diploid) intermingled with R. variabilis (allotetra- 
ploid). The synthetic di- and triploid hybrids generally clustered in 
an intermediate position between their parents, i.e., between the 
auricomus and the cassuhicus clusters. Although originating from 
different maternal progenitors, all synthetic F] hybrids were 
intermixed, forming a cluster of low internal heterogeneity (with 
the exception of two triploids (both from "H" population; 
Table 1), one of which grouped within the auricomus, another 
in the cassuhicus clade). In total, out of 501 nucleotide positions, 
34 were polymorphic in the directly sequenced ITS1-ITS2 
accessions (Figure 3). All 14 polymorphic sites detected in the 
ITS1-ITS2 of experimental hybrids were additive, either present- 
ing a combination of maternal [R. carpalicola or R. cassubicifolius) 
and paternal [R. notabilis) contributions (e.g., sites 52, 60 and 67), 
or identical to polymorphisms already existing in parental 
genomes (e.g., sites 26 and 97). The putative ancient hybrid 
Ranunculus variabilis displayed 15 polymorphic sites, three of 
which were shared with both di- and triploid synthetic hybrids. 
Interestingly, the three sexual species exhibited intraspecific 
polymorphisms, especially R. notabilis, which showed nine 
polymorphisms shared with R. variabilis. 

Inter- and intraindividual variability of ITS1-ITS2 inferred 
from cloning 

The average genetic identities within clone hbraries { — 
individuals) differed only slightly among the four cloned individ- 
uals and ranged from 0.998 (R. cassubicifolius 8472-1) to 0.996 {R. 
notabilis 7223-10). Averaged genetic similarities between the 
cloned individuals were highest between R. cassubicifolius and R. 
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The most frequent ITS variant of R. notabilis includes clones: 
NOB* and N06, N11-N13, N16, N18 



The most frequent ITS variant of R. cassubicifoliu s includes clones: 
S02* and S05-S16, S19, S20, X09, X12, X13, X16, X18, X19, X22 



Figure 4. NeighborNet analysis of cloned ITS1+ITS2 variants. Particular clones belong to: Ranunculus cassubicifolius ("S"-clones), R. notabilis 
("N"-clones), their synthetic hybrid ("X"-clones) and the putative hybrid, R. variabilis ("V"-clones). The spring leaf silhouettes illustrate the main 
phenotypic difference between the auricomus and the cassubicus morphotypes. Italic letters mark clones which exhibit non-compensatory base 
changes in the ITS1 or ITS2 secondary structures. Bootstrap values are given for the main clusters. 
doi:1 0.1 371 /journal.pone.01 03003.g004 



cassubicifolius x notabilis (0.999), and lowest between R. 
cassubicifolius and R. notabilis (0.973). Intraindividual ITSl- 
ITS2 variability and interindividual overlaps were visualized with 
NeighborNet (Figure 4). The splits diagram pointed out a high 
degree of intraindividual variability especially in the case of the 
apomict R. variabilis. The majority of R. variabilis clones 
clustered apart from a group of R. notabilis clones. However, 
several R. variabilis clones clustered in close vicinity to R. notabilis 
and even R. cassubicifolius x notabilis clones, suggesting the 
existence of similar ribotypes distributed over all three individuals. 
On the other hand, the R. cassubicifolius clones clustered apart 
from the R. notabilis/R. variabilis cluster. They formed one 
cluster together with the vast majority of the R. cassubicifolius x 
notabilis clones. As suggested from the splits graph (Figure 4), the 
most polymorphic individuals were the hybrids, i.e., R. variabilis 
and R. cassubicifolius x notabilis. Cloned ITS1-ITS2 amplicons 
from the four plant individuals exhibited 65 sites with nucleotide 



polymorphisms (i.e., 39 in ITSl and 26 in ITS2, Figure 5), of 
which 44 were not detected by direct sequencing. Within the 
entire ITS1-ITS2 region the synthetic and natural hybrid 
individuals shared three nucleotide polymorphisms (site 33, 431, 
.569, as marked with an asterisk in Figure 5). Furthermore, other 
shared nucleotide polymorphisms were detected for R. notabilis, 
its hybrid and R. variabilis (site 431) and for R. cassubicifolius, its 
hybrid and R. variabilis (site 569) (Figure 5). The ribotype 
network (Figure 6) points out a remarkable overlap among the 
4x R. cassubicifolius ribotypes ("S"-clones) and the ribotypes of its 
hybrid progeny R. cassubicifolius x notabilis ("X"-clones). 
Another prominent feature is the high amount of different ITS 1 - 
ITS2 ribotypes of R. variabilis ("V"-clones), where 95% of 
sequenced clones belonged to different ribotypes. None of these 
numerous ribotypes was completely identical to those from either 
R. notabilis ("N"-clones) or R. cassubicifolius x notabilis ("X"- 
clones). 
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Figure 5. Distribution of the ITS1+ITS2 nucleotide polymorphisms within the clone libraries. The polymorphic sites are represented by 
blue letters and marked with dots. Red letters marked with red asterisks represent polymorphic sites, which are shared by the synthetic hybrid of R. 
cassubicifolius and R. notabilis, their putative natural hybrid R. variabilis and any of the parental species. Double black dots mark polymorphisms which 
occur exclusively in single clones exhibiting non-compensatory base changes in the ITS1 or ITS2 secondary structures. 
doi:10.1371/journal.pone.0103003.g005 



Phenotypic variability within the "auricomus" 
morphotype 

Both di- and triploid synthetic hybrids (Figure 7a) resembled 
their (same) paternal parent, R. notabilis, having deeply dissected 
basal leaves [auricomus type). In contrast, both maternal parents 
(i.e., R. carpaticola, R. cassubicifolius) exhibited typical cassubicus 
morphology, characterized by non-dissected basal leaves (Fig- 
ure 2). The 2D landmark data (Figure 7b) revealed slight 
differences between R. notabilis, the synthetic hybrids and R. 
variabilis. The resulting scatter plot (Figure 7c) summarizes two 



most prominent morphological trends within the dataset of 95 
leaves. The main morphological trend (P' relative warp, RWl) 
describes variance at the leaf basis and accounts for 55.80% of the 
total morphological variability. This major morphological trend 
appears to reflect leaf size and does not separate the three groups. 
Only the second morphological trend (2""^ relative warp, RW2), 
describing a gradient from robust coarsely dissected to finely 
dissected forms, slightly differentiated R. variabilis and the 
synthetic hybrids, assigning R. notabilis an intermediate position 
in between. The mean phenetic distance (quantified as Procrustes 



PLOS ONE I www.plosone.org 



9 



July 2014 I Volume 9 | Issue 7 | e103003 



ITS Polymorphisms in Hybrids of Ranunculus 




°^"^^'//s 7223-10 (2X^ 

Figure 6. Ribotype network representing all cloned ITS1+iTS2 variants. Particular clones belong to: R. cassubicifolius ("S"-clones), R. notabilis 
("N"-clones), R. cassubicifolius x notabilis synthetic hybrid ("X"-clones) and R. variabilis ("V"-clones). The spring leaf silhouettes illustrate the main 
phenotypic difference between the auricomus and cassubicus morphotypes. Italic letters mark clones which exhibit non-compensatory base changes 
in the ITS1 or ITS2 secondary structures. 
doi:1 0.1 371 /journal.pone.01 03003.g006 

distance between two group centroids) was nearly the same R. variabilis and the synthetic hybrids accounted for nearly as 
between R. notabilis and R. variabilis (0.0536) as between R. twice as much (0. 1223). This was the sole case for which the mean 
notabilis and its synthetic hybrids (0.0738). The distance between morphological difference (based on the scores on 12 relative warps) 
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Figure 7. Geometric morphometric analysis of individuals exhibiting the characteristic "auricomu^' morphology. (A) An example of 
the sampled fresh plant material for DNA sequencing and geometric morphometric analyses. (B) 2D landmarks digitalization on the leaf outline. (C) 
Principal components analysis of the shape variables (i.e., Relative warps analysis, RWA) extracted from R, notabilis (blue squares), R. cassubicifolius x 
notabilis (red crosses) and R. variabilis (violet triangles) 2D landmark data. (D) Mean shapes of R. notabilis, R. cassubicifolius and R. variabilis 
reconstructed from each centroid (visualized as symbols surrounded by black-lined circles) of the three scatter clusters shovj/n in (C). The between- 
group differences were tested by permutation tests (lower left triangle, *** = p<0.001) and the distances between the mean shapes are expressed as 
Procrustes distances (upper right triangle). 
doi:10.1371/journal.pone.0103003.g007 



was statistically significant. Figure 7d visualizes the differences 
among the mean shapes of R. variabilis, R. notabilis and the 
experimental hybrids, illustrating the variation in the angle 
between the upper and lower part of the basal lateral lobes. 

ITS1-5.8S-IT52 polymorphisms in the secondary structure 
models 

The Ranunculus notabilis minimum energy ITS1-5.8S-ITS2 
secondary structure model (Figure SI) revealed 58% of the 241 
nucleotides involved in pairings in ITSl and 72% of the 212 
nucleotides in ITS2. In all directly sequenced individuals and 
clones, the conserved ITSl angiosperm motif 5'-GGCRY- 
(4- 7n)GYGYCAAGGAA-3' of [82] appeared as 5'- GGCGC- 
(GAUYG)GCGUCAAGGAA-3' in Ranunculus and contained 
no nucleotide substitutions. From the three 5.8S motifs that are 
also conserved across angiosperms, only one (M2) exhibited 
nucleotide substitutions (Figure 8), detected in three cloned 
amplicons. The conserved motif M3 (5'-TTTGAAYGCA-3') as 
described by [40], appeared as 5'-TTTGAACGCA-3'. In order to 
highlight and summarize the polymorphic sites detected by direct 
sequencing and cloning, they were mapped onto secondary 
structure models of 5.8S, ITSl and ITS2 rRNA transcripts 
reconstructed ior Ranunculus notabilis 5613-1 (Figure 8a, 9, 10). 
In the ITSl secondary structure we identified six hairpins (I- VI, 
Figure 9), in the ITS2 four hairpins (TIV, Figure 10). Those 
nucleotide substitutions which affect only one strand of the double- 
stranded helical regions are referred to as hemi-compensatory base 



changes (hCBC). Considering the directiy sequenced ITS1-ITS2 
accessions, out of the total number of 34 polymorphic sites, 12 
correspond to hemi-compensatory base changes and six base 
changes are not compensatory. The cloned amplicons yielded a 
total of 65 polymorphic sites in ITS1-ITS2. The additional DNA 
weblogos in Figures 9 and 10 illustrate all possible variants for 
each particular polymorphic site (as detected by direct sequenc- 
ing), and help to recognize non-compensatory (i.e., altering the 
secondary structure by missing pairing within any of the helical 
regions) vs. compensatory polymorphisms. Within the highly 
conserved 5.8S region, 13 sequences exhibited polymorphisms in 
single clones (Figure 8b) and might therefore represent putative 
non-functional ITS copies or cloning/PCR artifacts. 

Discussion 

Despite the tremendous popularity of the ITS nrDNA marker in 
angiosperm systematics, its variability on intraspecific or even on 
intraindividual level is still poorly known. The occurrence of ITS 
in multiple copies within a single genome requires more 
comprehensive sampling approaches to reconstruct reticulate 
phylogenies. Moreover, recovery of the whole range of individual 
variation can facilitate addressing questions regarding species 
evolutionary history. Here we contribute with the first detailed 
insight into the intraspecific ITS diversity of the Ranunculus 
auricomus species complex. Similar amount of recovered poly- 
morphisms across both ITSl (19%) and ITS2 (14%) regions was 
found in other polyploid complexes (e.g., [31,83]). 
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Figure 9. ITS1 secondary structure model for Ranunculus notabills 561 3-1 . Polymorphic sites detected in all studied taxa were highlighted in 
the secondary structure model for R. notabilis 5613-1. The DNA weblogo summarizes all ITS1 nucleotide polymorphisms detected by direct 
sequencing (red letters) and both direct sequencing and cloning (green letters). Bonds within the ITS1 secondary structure are shown for all positions 
which can be affected by either non-compensatory or compensatory (= hemi-CBCs) single-stranded nucleotide polymorphisms within helices. Italic 
letters mark non-compensatory base changes which are present only in a single clone. 
doi:1 0.1 371/journal.pone.01 03003.g009 



ITS polymorphism and morphology in the synthetic 
hybrids 

Our investigated reproduction system consists of a single 
apomictic microspecies, three sexuals with different ecology/ 
distribution areas and their synthetic hybrids. The diploid and 
triploid synthetic hybrids (Ti generation) were predominantiy 
sexual i^vitli a small j)(Tix'ntage of apomictic seeds in triploids [58]) 
and exhibited complete additivity in the ITS1-ITS2 polymor- 
phisms, reflecting genome contributions from their respective 
parents R. carpaticola Q, R. cassubicifolius Q and R. notabilis Cf. 
This pattern was already observed in the NeighborNet analysis of 
the directly sequenced individuals, where the hybrids occupied 
intermediate positions between their parental species, but is also 
evident in shared polymorphisms in the cloned data (Figure 5). 
When considering all cloning data, the vast majority of cloned 
ribotypes of synthetic hybrids were distinctiy shifted towards their 
maternal progenitor. This result is not surprising, as gametes of 
tetraploid R. cassubicifolius contributed two genome copies to the 
triploid offspring, whereas those of diploid -R. nolahilis contributed 
only one. Asymmetry in ITS copies in natural allopolyploidization 
events could have arisen from a higher number of molecular 
variation introduced by the tetraploid parent. Other authors 
[58,84] reported on a deviation from parental additivity after 
hybridization, the maternal dominance being attenuated after 
pol^ploidization as well. This uneven parental genome contribu- 
tion may be another reason why the putative natural hybrid, R. 
variabilis, does not match exactly the synthetic hybrids in the 
ribotype analysis, but appears to be more similar to R. notabilis. 
Nonetheless, the original hybridization event in nature could have 
happened with equal genome proportions between diploid R. 
cassubicifolius, which still occurs in the sympatric area with R. 
variabilis in one isolated population at the Austrian/Hungarian 
border [85]; Figure 1), and diploid R. notabilis. After hybrid 
origin, the early hybrid generations were probably sexual or only 
facultative apomictic, as in our experimental crosses [58]. 
Backcrossing of the first hybrid generations to the geographically 
closer parent, R. notabilis, or also concerted evolution of ITS 
copies in the first hybrid generations could have resulted in a 
higher genetic similarity of R. variabilis to R. notabilis. 

In contrast to other studies, referring to hybrids as either 
intermediate phenotypes or mosaics of their parents [3,46,86,87] 
our cloned ITS data underlined the remarkable incongruence 
between the genotype and leaf morphology of the Ranunculus 
cassubicifolius x nolahilis triploid hybrid. While including two 
genome copies from the maternal parent, the hybrid retained the 
clearly distinct paternal-like leaf morphology. Noack [88] 
described pentaploid hybrids that were more similar to the 
paternal progenitor rather than to the triploids which were 
morphologically intermediate. Further evidence against simple 
intermediacy was provided in [89] on Sorbns hybrid microspecies, 
observing leaves of mosaic, intermediate or parental-like features. 
As already discussed by [29,85], hybrid progeny in R. auricomus 
exhibits a continuous range of morphotypes ranging from extreme 
forms over intermediates to those completely identical with one of 
their parents. In the case of our synthetic hybrids, the observed 
parental-like leaf phenotype may also have an epigenetic basis. 
Results from our synthetic hybrids show that the parentage and 



evolutionary origin cannot be readily predicted from morpholog- 
ical data, while ITS polymorphisms reflect the genome additivity 

of the parental species. 

ITS polymorphism, reticulate evolution and mutation 
accumulation in the R. auricomus complex 

Compared to allozyme genetic identity within R. notabilis 
(0.965) andiJ. variabilis (0.940) assessed by [55], tiie ITS1-ITS2 
cDNA diflerentiation was in both cases considerably lower in these 
species (R. notabilis: 0.996, R. variabilis: 0.989). Although there 
was no particular ribotype shared between R. notabilis and R. 
variabilis, they clustered together. Liao et al. [90] recovered ten 
ITS polymorphisms in allotetraploid R. cantoniensis , a similar 
amount compared to our 13 polymorphisms detected in the 
allotetraploid R. variabilis. Intraindividual variability in ITS 
nrDNA can be maintained in apomictic hybrid species, as 
documented already by [30] in Amelanchier . Since apomictic 
reproduction generally results in seed formation without meiosis 
and syngamy, the homogenizing effects on ITS copies and 
concerted evolution are lacking. Hence, agamic lineages can 
preserve ITS polymorphisms from hybrid origin and can 
accumulate mutations independently in the polyploid genome. 
This so-called Meselson effect was in fact observed in hexaploid 
hybrids of the R. auricomus complex via large scale lUumina RNA 
sequencing [91]. Particularly in the case of allotetraploid R. 
variabilis, the observed amount of ITS polymorphisms may be 
indicative for lack of concerted evolution and post-origin mutation 
accumulation, a phenomenon frequently observed in apomictic 
plant polyploid species complexes [21,30]. 

Other than in the case of functional ITS copies, other authors 
[37] found most of the non-functional and chimeric ITS copies in 
a few million years old clade within Gnetum. We detected the 
largest number of putative non-functional copies (which contained 
polymorphisms in 5.8S region) in R. notabilis (20% of the clones) 
and in R. variabilis {16% of the clones). However, two 
polymorphisms in 5.8S were detected in the synthetic hybrid as 
well (10% of the clones). The ancient progenitor-derivative 
relationship between R. notabilis and R. variabilis, as proposed 
already by [5,5] based on allozyme data, was further supported by 
the observed overlaps in 10 shared polymorphic sites as inferred 
from direct sequencing. The high intraindividual ITS polymor- 
phism in R. notabilis individuals is unusual for a sexual species. 
Regarding the assessed split between R. notabilis and R. 
cassubicifolius about 0.9 Ma ago [54], such a high amount of 
ITS polymorphism is not expected. The most plausible explana- 
tion is that R. variabilis has introgressed into R. notabilis in early 
sexual hybrid generations, or also recentiy via backcrossing as the 
pollen donor, which is reasonable from the overlap of distribution 
areas (Figure 1) and the co-occurrence of the two taxa in some of 
the sites [55]. 

Secondary structure, compensatory base changes and 
non-functional ITS copies 

The amount of detected nucleotide pairings in RNA transcripts 
in the Ranunculus auricomus complex was higher compared to 
the observations by [65] on Asteraceae suggesting at least 20% of 
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ITS2 polymorphisms in the studied taxa detected 

direct sequencing - 4/30 

direct sequencing and cloning - 7/30 

cloning - 19/30 
ITS2 single-strand base changes within helices : 

compensatory - 421 , 428, 431 , 488, 514, 548, 

557, 571, 580, 599, 600 

non-compensatory - 455 ,519, 569, 576, 577, 

593, 601 
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Figure 1 0. ITS2 secondary structure model for Ranunculus notabllis 561 3-1 . Polymorphic sites detected in all studied taxa were highlighted 
in the secondary structure model for ft. notabilis 5613-1. The DNA weblogo summarizes all ITS2 nucleotide polymorphisms detected by direct 
sequencing (red letters) and both direct sequencing and cloning (green letters). Bonds within the ITS2 secondary structure are shown for all positions 
which can be affected by either non-compensatory or compensatory (= hemi-CBCs) single-stranded nucleotide polymorphisms within helices. Italic 
letters mark non-compensatory base changes which are present only in a single clone. 
doi:l 0.1 371/journal.pone.01 03003.g01 0 



ITSl and 38% of ITS2. N(;v(:rth(:l(\ss, iu'ith(^r among the sexual 
species nor between them and the apomict R. variabilis did we 
recover any CBCs. The importance of these structural features for 
species concepts is still under debate [33,34]. The strongly reduced 
fertility of our synthetic hybrids [92] suggests that their sexual 
parents are distinct biological species. Hence, the relevance of 
CBCs for delimitation of species, and for reconstructing reticulate 
evolution and interspecific hybridization remains unclear. How- 
ever, the occurrence of the numerous compensatory single-strand 
nucleotide polymorphisms (i.e., hemi-CBCs) might be a hint for an 
ongoing process resulting in a complete CBC during further 
divergence of taxa. More than a half of the hemi-compensatory 
base changes in ITS 1 -ITS2 were present exclusively in the cloned 
ampKcons. A suspected artificial nature can be excluded because 
the observed substitutions retain the nucleotide bonds and the 
secondary structures. In contrast, the majority of non-compensa- 
tory substitutions were present in single clones only and hence, 
may represent PCR artifacts. Furthermore, the secondary 
structure analysis of the 5.8S region indicated no substitutions in 
the directly sequenced amplicons, but numerous substitutions in 
single clones. Sixteen percent of the analyzed clones showed 
polymorphisms in the 5.8S, belonging to 15 different variants, an 
observation close to [41] from Mammillaria. In contrast to the 
mutations present in the DNA template, the nucleotide substitu- 
tions in otherwise conser\'ed 5.8S might have been introduced 
during PCR reaction and E. coli cell division [93-95] . Neverthe- 
less, Zaveska and collaborators [21] identified in Taraxacum 
considerable polymorphisms in 5.8S (3.7%), considering all these 
ITS copies to be functional with regard to their secondary 
structures. Similarly, we neither observed structural changes 
within 5.8S introduced by the substitutions nor did we recover 
any recombinant/chimeric sequences through recombination 
detection algorithms or via NeighborNet analysis. Additionally, 
the patterns of methylation-induced substitutions or GC content 
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